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The nonlinear pulse propagation in an optical fibers with varying parameters is investigated. The 
capture of moving in the frequency domain femtosecond colored soliton by a dispersive trap formed in 
an amplifying fiber makes it possible to accumulate an additional energy and to reduce significantly 
the soliton pulse duration. Nonlinear dynamics of the chirped soliton pulses in the dispersion 
managed systems is also investigated. The methodology developed does provide a systematic way 
to generate infinite "ocean" of the chirped soliton solutions of the nonlinear Schrodinger equation 
(NSE) with varying coefficients. 
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I. INTRODUCTION 

In 1973 Hasegawa and Tappert [|] showed theoretically that an optical pulse in a dielectric fibers forms an envelope 
solitons, and in 1980 MoUenauer, Stolen and Gordon j^] demonstrated the effect experimentally. This discovery 
is significant in its application to optical communications. Today the optical soliton is regarded as an important 
alternative for the next generation of high speed telecommunication systems. 

The theory of NSE solitons was developed for the first time in 1971 by Zakharov and Shabad |^ . The concept of the 
soliton involves a large number of interesting problems in applied mathematics since it is an exact analytical solution 
of a nonlinear partial differential equations. The theory of optical solitons described by the nonlinear Schrodinger 
equation has produced perfect agreement between theory and experiment Q . 

In this paper we present mathematical description of solitary waves propagation in a nonlinear dispersive medium 
with varying parameters. 

The soliton spectral tunneling effect was theoretically predicted in 1^ . This is characterized in the spectral domain 
by the passage of a femtosecond soliton through a potential barrier-like spectral inhomogeneity of the group velocity 
dispersion (GVD), including the forbidden band of a positive GVD. It is interesting to draw an analogy with quantum 
mechanics where the solitons are considered to exhibit particle-like behavior. The soliton spectral tunneling effect 
^>~,' also can be considered as an example of the dynamic dispersion soliton management technique. In the first part 
of the paper we will concentrate on the problem of femtosecond solitons amplification. We will show that spectral 
Q-i inhomogeneity of GVD allows one to capture a soliton in a sort of spectral trap and to accumulate an additional 
^ , energy during the process of the soliton amplification. In the second part we will consider the problem of the short 
soliton pulse propagation in the nonlinear fiber with static non- uniform inhomogeneity of GVD. The methodology 
developed does provide a systematic way to generate infinite "ocean" of the chirped soliton solutions of NSE model 
^ I with varying coefficients. 
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II. FEMTOSECOND SOLITON AMPLIFICATION 



It is well known that due to the Raman self-scattering effect [|j (called soliton self-frequency shift the central 
femtosecond soliton frequency shifts to the red spectral region and so-called colored solitons are generated. This effect 
decreases significantly the efficiency of resonant amplification of femtosecond solitons. The mathematical model we 
consider based on the modified NSE including the effects of molecular vibrations and soliton amplification processes 
(see details in Q): 
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^,^^ + 2f,S^ + Q^\tP\', and, + P(l + *7aAr!) = z^-, (2) 

As numerical experiments showed the GVD inhomogeneity as a potential well allows one to capture a soliton in a 
sort of spectral trap. Figure |^ shows the nonlinear dynamics of the soliton spectral trapping effect in the spectral 
domain. As soliton approaches the well, it does not slow down but speeds up, and then, after it has got into the well, 
the soliton is trapped. There exists a long time of soliton trapping in internal region of the well .This effect opens 
a controlled possibility to increase the energy of a soliton. As follows from our computer simulations the capture of 
moving in the frequency space femtosecond colored soliton by a dispersive trap formed in an amplifying optical fiber 
makes it possible to accumulate an additional energy in the soliton dispersive trap and to reduce significantly the 
soliton pulse duration. 



III. DISPERSION MANAGEMENT: CHIRPED SOLITONS 



Let us consider the propagation of a nonlinear pulse in the anomalous (or normal) group velocity dispersion fiber 
of length Zi. The complex amplitude q of the light wave in a fiber with variable parameters D2{Z) , N2{Z) and r(Z) 
is described by the nonlinear Schrodinger equation 



(3) 



Theorem 1. Consider the NSE (3) with varying dispersion, nonlinearity and gain. Suppose that Wronskian 
W[N2,D2] of the functions N2(Z) and D2(Z) is nonvanishing, thus two functions N2(Z) and D2(Z) are linearly inde- 
pendent. There are then infinite number of solutions of Eq. (3) in the form of Eq.4 



q{Z,T) = 



lD2{Z) 
N2{Z) 



P{Z) Q [P{Z) ■ T] cxp 



K{Z )dZ 



(4) 



where function Q describes fundamental functional form of bright Q — sech(P(Z)r) or dark Q = th{P{Z)T) NSE 
solitons and the real functions P(Z), D2{Z), N2{Z) and T{Z) are determined by the following nonlinear system of 
equations : 



1 dPjZ) 
P^{Z) dZ 



D2iZ)^0 



^D2{Z)P{Z) 



W[N2{Z),D2{Z)] 
2D2{Z)N2{Z) 



T{Z) 



(5) 



Theorem 2. Consider the NSE (3) with varying dispersion, nonlinearity and gain. Suppose that Wronskian 
W[N2,D2] of the functions N2(Z) and D2(Z) is vanishing, thus two functions N2(Z) and D2(Z) are linearly dependent. 
There are then infinite number of solutions of Eq. (3) of the following form Eq. 6 



q{Z, T)^C P{Z) Q [P{Z) ■ T] cxp 



K{Z )dZ 



(6) 



where function Q describes the fundamental form of bright (or dark) NSE soliton and the real functions P(Z), D2(2'), 
'H2{Z) and r(Z) are determined by the following nonlinear system of equations : 



D2{Z) 



1 dP{Z) 
P^'{Z) dZ 



1 1 dP{Z) 
2P dZ 



N2{Z) = D2{Z)/C'' 



(7) 



The function P(Z) is required only to be once-differentiable, but otherwise arbitrary function, there is no restrictions. 

To prove Theorems 1 and 2 we first construct a stationary localized solution of Eq. (3) by introducing Kumar- 
Hasegawa's quasi-soliton concept [|| Q through 
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ID2{Z) 
N2{Z) 



P{Z) Q [P{Z) ■ T] cxp 



K{Z )dZ 



(8) 



where 1)2(2), N2{Z), P{Z) and K{Z) are the real functions of Z. Substituting expression (8) into Eq. (3) and 
separating real and imaginary parts we obtain the system of two equations 



E-^-n\Z)] Q = 



(9) 



Where 

S{Z, T) = P{Z)T ; II = T— ; — ^ P{Z) (11) 

In Eq. (9) the parameters E and O are 'the energy' and 'frequency' of ordinary quantum mechanical harmonic 
ocsillator 

^\Z) = (t^Z) S + ^^^^0 ' ^^^^ ^ -KiZ)/P^iZ)/D2iZ) (12) 

Eq. (9) represents the nonlinear Schrodinger equation for the harmonic ocsillator. As must be in the case of Hamil- 
tonian system Eq. (9) may be written in the form 



H = 



-sign(i:>2) 



dQ 
dX 



E-^-V?{Z) 



dX 



(14) 



The derivative in (13) is functional derivative. For the first time this equation was solved numerically by Kumar and 
Hasegawa in and gave rise a new concept of quasi-solitons ||l^;|ll| • Now we make the important assumption about 
the solution of Eq. (9). 

Let us consider the complete nonlinear regime when Eq. (9) represents the ideal NLS equation, i.e. we will allow 
n{Z) = , then from (12) follows that 



1 dPiZ) 
P^{Z) dZ 



+ D2{Z) = Q 



(15) 



We now look for a solution of Eq. (10) which satisfies the condition (15). Substituting the expression (15) and 
relations (11) into Eq. (10) we obtain 



^D(Z)P(Z)+^ ^ 1 1 ^^-rm 

- -D2{Z)P[Z) + 2^^(^) Qz ' 2N2{Z) dZ - ^^^^ 



(16) 



Using notation 



W{N2,D2} 



dD2{Z) dN2{Z) 

J>'2 — U2- 



dZ 



dZ 



one can obtain the soliton solution of Eq. 3 in the form of the chirped solitons Eqs. 4-5 and Eqs. 6-7.. Consequently, 
we have found the infinite "ocean" of solutions. The methodology developed does provide a systematic way of new 
and new chirped soliton solutions generation. 
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IV. DIFFERENT REGIMES OF SOLITON MANAGEMENT 



Lemma 1: Soliton GVD management. Consider the NSE (3) with constant nonUnear coefScient N2 = const and 
with varying along Z-coordinatc GVD parameter. Suppose that dispersion management function is known arbitrary 
analytical function :D2(Z)=$(Z) . The function ^{Z) is required only to be once-differentiable and once integrable, 
but otherwise arbitrary function, there is no restrictions. There are then infinite number of solutions of Eq. (3) of 
the form of the chirped dispersion managed dark and bright solitons Eq. 4, where the main functions P and F are 
given by 

i..(Z)=$(Z);P(Z) = -^^-^.l^ (17) 

1 ^jZ) 1 1 d<i>iZ) 

'^^> 2[C-mZ)dZ] 2<^{Z) dZ ^ ' 

Lemma 2: Soliton intensity management. Consider the NSE (3) with constant nonlinear coefficient N2 = const and 
with varying along Z-coordinate the dispersion and gain. Suppose that intensity of the soliton pulse is determined by 
the known management function: D2(Z)P^(Z)=0(Z), where the function Q{Z) is required only to be once-differentiable 
and once integrable, but otherwise arbitrary function, there is no restrictions. There are then infinite number of 
solutions of Eq. (3) of the form of the chirped dispersion managed dark and bright solitons Eq. 4 with parameters 
given by 



D2{Z)P\Z) = e{Z) ; P{Z) = - j e{Z)dZ + C ; D2iZ) 



Q{Z) 



[C-jQ{Z)dZy 



2[c- jeiz)dz] 2e{z) dz ^ ' 

Lemma 3: Soliton pulse duration management: optimal soliton compression. Consider the NSE (3) with constant 
nonlinear coefficient N2 = const and with varying along Z-coordinate the dispersion and gain coefficients. Suppose 
that pulse duration of a soliton is determined by the known analytical fimction: P(Z)=T(Z), where the fimction T{Z) 
is required only to be two-differentiable , but otherwise arbitrary function, there is no restrictions. There are then 
infinite number of solutions of Eq. (3) of the form of the chirped dispersion managed dark and bright solitons Eq. 4 
with the main parameters given by 

Lemma 4: Soliton amplification management: optimal soliton compression. Consider the NSE (3) with constant 
nonlinear coefficient N2 = const and with varying along Z-coordinate the dispersion and gain coefficients. Suppose 
that the gain coefficient is determined by the known control function: T{Z)=A{Z), where the function A{Z) is required 
only to be once integrable . but otherwise arbitrary function, there is no restrictions. There are then infinite number 
of solutions of Eq. (3) of the form of the chirped dispersion managed dark and bright solitons of the Eq. 4 where 



|P(Z)|=exp 



2A{Z )dZ dZ 



(21) 



IP .^z. exp(/2A(Z)ciZ) 

' " exp [/exp {j2K{Z")dZ") dZ'] ^ ' 

Lemma 5: Combined dispersion and nonlinear soliton management. Consider the NSE (3) with varying nonlinear 
coefficient N2 {Z) and with varying along Z-coordinate the dispersion and gain coefficients too. Suppose that Wronskian 
W[N2,D2] is vanishing, or that the functions N2(Z) and D2(Z) are linearly dependent. Suppose also that the function 
D2(Z) is determined by the initial control fimction D2(Z)=S(Z), where the fimction S(Z) is required only to be once 
integrable, but otherwise arbitrary function, there is no restrictions. There are then infinite number of solutions of 
Eq. (3) of the form of the chirped dispersion managed dark and bright solitons of the Eq. 6 where 
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P{Z) = 



[C-jEiZ)dZ] 



; N2{Z) = D2{Z)/C^ 



(23) 



1 



E{Z) 



2[C-jE{Z)dZ] 



(24) 



The analytical solutions for the different regimes of the main soliton parameters management (intensity, pulse duration, 
amplification or absorption ) in the case of W[N2,D2]=0 can be obtained by using theorem 2. 

Let us consider some examples. The case of r(Z) = and N2(Z)=N2(0) corresponds to the problem of ideal GVD 
soliton management. The soliton solution in this case is: 



(?(Z,T) = -TjN^ ^^^(0)exp( jZ ) sech [T?Texp(CZ)] 



C 



exp 



C 1 
^iT^— exp(CZ) -i-ri^Ze^p{CZ) 



(25) 
(26) 



q{Z,T) = r^N^ '/'(0)exp( |z ) th [r;rexp(CZ] 
C 

exp iT^—exp{CZ)-ir]^Zei<p{CZ) 



(27) 
(28) 



Here T and Z are ordinary variables and C is arbitrary constant. If we use the expressions D2(Z)=constant and 
N2=N2(0) then we obtain the following solutions of Eq. (3) in the form of hyperbolically growing ideal bright and 
dark solitons (for the first time reported in p^ , p^ 



,(Z,r)^-^^^";(»\ sech 

' ' (i-2r(o)z) 



xT 



(1 -2r(o)z) 



exp 



r2r(o) 



(i-2r(o)z) 2(i-2r(o)z) 



(29) 



,(z,T)^,^^!r?) th 



(1 - 2r(o)z) 



exp 



r2r(o) 



x'z 



[1 - 2r(o)z) 

In the case of T{Z) = Gq and N2=A''2(0) the solution of Eq. 3 is given by: 

Q{P{Z)T) = riN:^^''^{Q) sech [riP{Z)T] 

Q{P{Z)T) = r,N-"\Q) th [vP{Z)T\ 



P{Z) = -P(0)exp(— (exp(2GoZ) - 1)) 



(i-2r(o)z) (i-2r(o)z) 



(30) 

(31) 
(32) 

(33) 



D2{Z) - i?2(0) exp(2Go^ - — (exp(2Go^) - 1)) 



When GVD is a hyperbolically decreasing function of Z 



(34) 



D2{Z) 



1 + /3Z 

then from Lemma 1 follows the explicit soliton solution in the form of Eq. 4 

1 



iln(l + /3Z) 



(35) 



(36) 
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T{Z) 



1 



2(1 + /3Z) 



1 - ln(l + I3Z) 
1 - i ln(l + (iZ) 



(37) 



Let us consider the soliton intensity management problem. Chirped soUton pulse of Eq. 3 with the constant 
intensity can be obtained by using Lemma 2 



P{Z) ^ -CZ - 1; D2{Z) = C/(l + CZf- T{Z) = -C/2/(l + CZ) 



(38) 



Let us consider some periodical chirped soliton solutions of Eq. 3. Suppose that the soliton intensity varies periodically 
as 



D2{Z)P'^{Z) = Q{Z) = 1 + (5sin2"Z 
Then soliton solution in the case of n=2 is determined by Eq. 4 with parameters: 

' iZ sm2Z sin4Z 



D2{Z) = Q{Z)IP'{Z); P{Z) = C 



Z + 5 



32 



(39) 



(40) 



1 {l + 5s\n^Z) 1 2 sin 2Z sin^ Z 

2C-[Z + <5(f -2H^ + 2i|M)] +2 (l + ^sin^Z) ^ ' 

Let us consider some periodical solutions of Eq. 3 in the case of the linearly dependent parameters of the media. The 
simplest solution of Eq. 3 in the form of Eq. 6 is: 

P{Z) = T{Z) = - (1 + <5sin2 Z) ; N2{Z) = D2{Z) = '^'^"^f ; (42) 

(l + (5sin^ Z) 

, , (5 sin2Z 

2 [^l + d sm Z) 

The next periodical soliton solution is given by 

D2iZ) . N2iZ) = cosZ; P(Z) = ^j^^^y r(Z) = ^j^^^ (44) 

The main soliton features of the solutions given by theorem 1 and theorem 2 were investigated by using direct 
computer simulations. We have investigated the interaction dynamics of particle-like solutions obtained, their soliton- 
like character was calculated with the accuracy as high as 10~^. We also have investigated the influence of high-order 
effects on the dynamics of dispersion and amplification management. As follows from numerical investigations elastic 
character of chirped solitons interacting does not depend on a number of interacting solitons and their phases. Figure^ 
shows the computer simulation dynamics of three hyperbolically growing solitons Eq. 29. NSE solution with periodic 
dispersion coefficient is shown in Figure ^. Here the dispersion management function is 

D2iZ) = l + Ssm^{Z) (45) 

and the soliton solution is given by Eqs. 17-18. In Figure ^ parameters C=200 and S = —0.9. Figure ^ represents the 
two dispersion managed solitons interaction in the case of equal phases and in Figure |^ the interaction dynamics of 
two solitons is shown in the case of opposite phases. Figure ^ shows the intensity managed solitons dynamics of the 
form presented by Eq. 38. Figures show the nonlinear propagation and interaction of the dispersion and nonlinear 
managed solitons of Eqs. 42-43. The main parameters in computer simulations were C=200; S = ±0.9. Figure |l^ 
illustrates the dynamics of the fission of the bound states of two hyperbolically growing solitons Eqs. 29-30 produced 
by self-induced Raman scattering effect given by Eqs. 2-3. This remarkable fact also emphasize the full soliton features 
of solutions discussed. They not only interact elastically but they can form the bound states and these bound states 
split under perturbations. The possibility to find the plethora of soliton solutions in the case of strong dispersion 
management is reported in the recent paper of Zakharov and Manakov llj]. 
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FIG. 1. Femtosecond soliton spectral trapping effect. 

FIG. 2. Mutual interaction of three hyperbolically growing chirped solitons of Eq. 29 in the case of equal amplitudes and 
phases. 

FIG. 3. Evolution of the chirped dispersion managed solitary wave of Eqs. 17 and 18 as a function of the propagation 
distance. Dispersion managed function is periodic of the form Eq. 45. Input conditions : C=200 and 5=-0.9. 

FIG. 4. Two dispersion managed solitons of Eqs. 17-18 and 45 interaction for the case of equal phases. Input conditions: 
C=200 and S=-0.9. 



FIG. 5. Two dispersion managed solitons of Eqs. 17-18 and 45 interaction for the case of equal phases. Input conditions: 
C=200 and (5=-0.9. 



FIG. 6. Two intensity managed solitons of Eq. 38 interaction for the case of zero initial group velocities. 



FIG. 7. Evolution of the chirped solitaty wave of Eqs. 39-41 as a function of the propagation distance. Input conditions: 
6=0.9 and group velocity Vb=10. 



FIG. 8. Evolution of the chirped solitary wave of Eqs. 39-41 for the case: S =-0.8 and group velocity Vb=2.0. 

FIG. 9. Soliton dispersion trapping effect in the presence of the linearly dependence between the nonlinearity and GVD 
parameters. 

FIG. 10. Decay of high-order hyperbolically growing solitons in the presence of third-order dispersion and Raman 
self-scattering effects. 
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